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Abstract 

Generalized Wess-Zumino models which admit topologically non-trivial 
BPS saturated configurations along one compact, spatial dimension are inves- 
tigated in various dimensions of space-time. We show that, in a representative 
model and for sufficiently large circumference, there are BPS configurations 
along the compact dimension containing an arbitrary number of equidistant, 
well-separated domain walls. We analyze the spectrum of the bosonic and 
fermionic light and massless modes that are localized on these walls. The 
masses of the light modes are exponentially suppressed by the ratio of the 
distance between the walls and their width. States that are initially local- 
ized on one wall oscillate in time between all the walls. In (2+1) dimensions 
the "chirality" of localized, massless fermions is determined. In the (1+1)- 
dimensional case we show how the mass of certain classically BPS saturated 
solitons is lifted above the BPS bound by instanton tunneling. 
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Introduction 



Recently, there have been various proposals entertaining the idea that our d = 3 + 1 
dimensional world is located on a brane in a higher dimensional space-time. In 
such scenarios, the standard model fields are confined to the brane, whereas gravity 
permeates the bulk. The idea of compactification by localization of zero modes on 
defects of a higher dimensional field theory was first advocated in Ref. [|1|. 

In Refs.p], H, H, the hierarchy problem is addressed through the notion of n 
compact (large) extra dimensions of size R < 1 mm. In this scenario the Planck scale 
M4 of our world is not fundamental but can be derived from a higher dimensional 
scale M 4+n of the order of the electroweak scale My/ by means of the following 
power-law conversion: 

M 2 A ~ MlX^R n . (1) 

However, there is still a large hierarchy of the fundamental scales M4 +n , 1/R. 

In contrast, the authors of Ref. || suggest to solve the hierarchy problem in a 
rather different way. The set-up is a five-dimensional universe consisting of one three- 
brane constituting our world and another, hidden three-brane, which are separated 
by a (small) distance nr c . Assuming a special, non-factorizable metric for the five 
dimensional space-time, there is practically no conversion between (fundamental) 
Planck scale M 5 and the Planck scale M 4 of our world. Scales M belonging to the 
electroweak regime are given in terms of M 5 as 

M = e~ 2nkrc M 5 , (2) 

where the fundamental scales k,l/r c have no large hierarchy (kr c ~ 50). 

In addition, in Ref. |7]] it was proposed that the fermion hierarchy in the standard 
model may be generated by localizing different families on spatially separated walls 
along a compact dimension. Thereby, the localized light modes obtain their masses 
by coupling to a Higgs condensate which falls off exponentially away from the wall 
on which the Higgs field is localized. 

The purpose of this paper is to study the idea of "brane worlds" in a very simple 
super symmetric context, namely Wess-Zumino models. One and the same field is 
used to create domain walls (which play the role of the branes) and the fermionic 
and bosonic light modes localized on them. The latter play the role of matter in the 
lower dimensional universes on the walls. We will assume the geometry of space- 
time as given, and we do not investigate the issue of radius stabilization, although 
the model we discuss in fact provides a realization of "brane-lattice-crystallization", 
proposed in Ref.|5j as a solution to this problem. The models we consider feature 
several distinct topological sectors. We do not address the question why a particular 
sector is selected. Obviously, such a limited approach poses severe restrictions: 
there is no gravity, no gauge fields, and we have to work in lower dimensions. The 
advantages are transparency, the ability to look at different issues in isolation, and 
the potential to perform explicit calculations. We will be able to address issues such 
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as spontaneous supersymmetry breaking, the chirality of the localized fermionic 
modes, and the spectrum of the localized light modes. 

For concreteness, we choose to work with a specific model that has the property 
that we desire; static configurations of well separated domain walls in a compact 
space dimension. We stress that it is just a toy model; the conclusions we draw 
should not, and do not, depend on its specific details. 



The model 

To be specific, we consider a generalized Wess-Zumino model in d — 3+1 dimensions 
with one chiral superfield 0.Q The superpotential of this model is given by 

W = In (0 + 1)+ In (0-1) , (3) 

and the Kahler potential takes the canonical form. We also consider this model 
reduced to d = 2 + 1 and d = 1 + 1 dimensions. 

Since the coordinate z is compact, the field obeys periodic boundary conditions 
4>(x) = 0(x + L), where L is the circumference of the compact dimension. In general, 
static, stable solutions to the second-order equations of motion may or may not be 
BPS saturated. If a BPS saturated solution is selected as the vacuum of the model, 
only 1/2 of the SUSY generators is spontaneously broken. We will focus mostly on 
the BPS saturated solutions, but static, stable, non-BPS solutions also exist, and 
we will comment on them in passing. 

Since the presence of a non-vanishing (1, 0) central charge Z in the SUSY algebra 
is necessary for the existence of BPS domain walls, it was argued in Ref. ]14| that, 



in generalized Wess-Zumino models with a compact dimension, the superpotential 
must be a multi-branch function and that the manifold M on which the fields are 
defined must admit non-contractable cycles. The model with the superpotential W 
of Eq. (P) fulfills these criteria. The manifold M in this case is a plane with two 
punctures a the branch points = ±1. 

As we also have occasion to consider the model in less than d = 3 + 1 dimensions, 
and as the scalar components of the M = 1 super multiplets in lower dimensions are 
real, it is convenient to decompose the field according to = l/\/2(0i + 202) with 
real fields 0i and 02- By definition, the form of the scalar potential is not altered 
by the procedure of dimensional reduction. The BPS equation and its solutions are 
the same in each number of space-time dimensions that we consider. Therefore, 
we will first characterize all periodic solutions to the BPS equation; armed with 
this knowledge, we will subsequently investigate various implications in different 
dimensions. Of course, a solution along a compact dimension is constrained by the 
requirement that its period A fits an integer N times on the circumference L. 



: We denote the superfield and its lowest component by the same symbol 
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The scalar potential of the model is given by 



_ dW dW 



+ 



+ 1 0-1 



(4) 



This potential is invariant under the transformations 0i — > — 0i and 02 - > — 02- It 
has a supersymmetric vacuum, V = at = 0, and also a run-away vacuum at 
|0| — ► oo. In the vacuum at the origin, the mass of the field is m = 2. In addition, 
there are two saddle points at = ±2, and two poles at = ±1. As we will discuss 
in detail below, the solutions to the BPS equation either wind around one of these 
poles, or around both. 

At this point we briefly pause to discuss radiative corrections to the classical 
BPS solutions and the spectrum of the localized modes. Obviously, the model with 
the superpotential Eq. (|3[) is not renormalizable in d — 3 + 1 dimensions. Clearly, it 
should be regarded as an effective theory, describing the low energy dynamics of a 
fundamental theory below an energy scale M. The field is supposedly the only 
relevant degree of freedom below this scale. 

In fact, in order for an effective theory to be useful, it should be in the weakly 
coupled regime. To that end, the superpotential in Eq.(||) can be slightly modified 
by the introduction of a coupling constant g, 

^ = ln(s0+l) + ln(s0-l) , (5) 

so that for g<l the model is certainly weakly coupled. 

The tension of the domain walls is T = 2tc, independent of g. The scale T 
should be much smaller than the scale M for consistency of the effective theory 
approach. Moreover, the mass in the supersymmetric vacuum at the origin now 
becomes m = 2g 2 , and the mass scale m pzm of the pseudo-zero modes localized on 
the wall is exponentially suppressed by the ratio of the distance between the walls 
and their width. We thus find the folowing hierarchy of scales 

M > T > m > m pzm . (6) 

It is the underlying fundamental theory that must create this hierarchy. As a conse- 
quence, it is possible to describe the dynamics in different energy ranges by various 
effective theories. 

Radiative corrections to local quantities such as the energy density profile of 
the domain wall can therefore be calculated in perturbation theory as described in 



Refs.[T7, 18, 13]. The scale M forms a physical ultra-violet cut-off in this calculation, 
and the size of the corrections is controlled by g. 

Moreover, the dynamics of the zero modes and pseudo zero modes that are 
localized on the walls can be described by a weakly coupled, d = 2 + 1 dimensional, 
effective theory. In this case, the scale m forms the physical cut-off. Above this 
scale, more massive modes become accessible and should be incorporated in the 
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theory. Radiative corrections can be calculated in a perturbative expansion and are 
controlled by powers of the coupling constant g. 

As radiative corrections are not the issue we investigate in this paper, we will 
just take g — 1 in what follows. Our results can be trivially generalized for other 
values of g. In addition, we will not speculate about the nature of an underlying 
fundamental theory that could give rise to the superpotential in Eq.(|]). 

The paper can be outlined as follows: In Section [I] we study and classify all 
periodic solutions to the BPS equation. We show that there are BPS saturated 
configurations that contain well separated domain walls. In Section § we investi- 
gate such configurations with a sequence of equidistant walls in (3 + l)-dimensions. 
Apart from the bosonic and fermionic zero modes associated with the spontaneous 
breaking of translational invariance and half of supersymmetry, there are other light 
modes with exponentially suppressed masses. The sequence of domain walls can 
be considered as a crystal, and the spectrum of the pseudo-zero modes is analyzed 
using tight binding methods. We then study to what extend the domain walls can 
be considered as parallel universes. After integrating out the heavy modes, the 
(2 + l)-dimensional theory of the light modes has M = 1 supersymmetry. In the 
multi-wall background, light states that are localized on one particular wall are no 
longer mass eigenstates. This leads to oscillation phenomena. Section [3] investigates 
the same model in (2 + 1) dimensions. The interesting new feature is that we are 
now able to study the "chirality" (whether they are left or right moving) of the 
fermionic zero modes localized on the domain walls (there are no chiral fermions in 
d — 2 + 1 dimensions). We show that in the BPS saturated background there are 
always one left moving and one right moving fermion. We argue that this continues 
to be the case even if the M = 2 supersymmetry is explicitly broken to M = 1, 
in which case only one SUSY generator is spontaneously broken. In Section |] we 
study the model in d — 1 + 1 dimensions. What is different here is that the domain 
walls with infinite mass in the higher dimensions now become solitons with a finite 
mass. We show that in particular topological sectors, and for sufficiently large cir- 
cumference of the compact dimension, the classical ground state is BPS saturated, 
but non-perturbatively the ground state energy is lifted above the BPS bound. This 
phenomenon was first observed in Ref. |2U|. The ground states are connected by 



instanton tunneling. The tunneling probability and the energy shift of the ground 
state are obtained as a function of the circumference of the compact dimension. The 
Appendix contains technicalities concerning the fermionic zero modes. 

1 Periodic solutions to the BPS equation 

In this section we study periodic solution to the BPS equation, 



dx 
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where the phase e takes the value ±z for periodic solutions. The two possible 
choices for the phase correspond to whether the solutions wrap around the poles 
clockwise or counter clockwise in the complex plane. Solutions for one choice can 
be obtained from solutions with the opposite choice of the sign by the transformation 
x — > —x, mapping walls into anti-walls and vice versa. For definiteness, we will from 
here on consider the BPS equation with e lS = +i. It is well known that solutions of 
the BPS equation have a "constant of the motion" / = Imfe -4 ' 5 ^]. This constant of 
the motion can in principle be used to mark the individual solutions. Here, instead 
of J, we use the related quantity 



to label the solutions, which takes a simpler form. The "constant of the motion" k 
is positive semi-definite. It only vanishes at the poles <f> — ±1. Therefore, there are 
solutions for each value k > 0. Since the BPS equation is a first order differential 
equation, solutions, when viewed as trajectories in the complex plane, do not 
intersect each other or themselves, except possibly at the supersymmetric vacuum 

= 0. 

The target space is a plane with two punctures at the location of the poles of 
the scalar potential. The trajectories of solutions to the BPS equation fall into 
three homotopy classes of non-contractable cycles in the complex plane.0 For any 
given value of k in the range < k < 1 there are two distinct solutions to the BPS 
equation. These two solutions are mapped onto each other by the transformation 

01 — ► — 0i) which is a symmetry of the scalar potential. One solutions winds around 
the pole at = 1, whereas the other winds around the pole at = —1. We will 
denote these two homotopy classes as Class IA and Class IB, respectively. Solutions 
in a third homotopy class, which we refer to as Class II, are obtained for k > 1. For 
each value of k > 1 there is only one solution, which winds around both poles of the 
scalar potential and is invariant under the transformation 0i — > — 0i. 

There is a critical solution for k = k c = 1 which separates solutions in the three 
homotopy classes in the complex plane. This critical solution contains the point 
= 0, where the scalar potential has a minimum. The wavelength of this solution 
is therefore infinite. It represents a domain wall that interpolates between the same 
supersymmetric vacuum at x — > — oo and x — > +oo, winding either around the pole 
at = —1 or the pole at = +1. We will refer to these walls as "left" and "right" 
walls respectively. Examples of solutions in each of the homotopy classes are shown 
in Fig. (|TJ) . At k = 4 the cycle passes through the two saddle points at ± 1. For 
the cycles in Class II, the wavelength tends to infinity both in the limit k j 1 and 
k — > oo. In between, the wavelength reaches a minimum which we numerically 
determined to take the value A m = 4.96 for k m = 1.42. 

2 Here, we are actually classifying the functions /fc(s) = 0fc(sA(fc)) (s e [0, 1]), where A(fc) is the 
wavelength of the solution labeled by k. 




(8) 
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Figure 1: BPS solutions in the complex plane. The critical solution is denoted 
by the dashed line. The dotted line indicates the Class II solution with minimum 
wavelength. Inside the critical solution, a Class IA and a Class IB solution wind 
around their respective poles. Also shown is a Class II solution passing through 
= 1, where the scalar potential has a saddle point. Two other Class II solutions, 
one with a value of k near one (encompassing the critical solution), the other with 
a larger value of k (outermost cycle) are indicated. 

For very small values of k, k <C 1, the cycles are small circles centered around 
the poles in the complex plane. The wavelength of such solutions is approximately 
A = nk/2, and the energy density is approximately constant. For large k, k 1, 
the cycles are approximately large circles centered around the origin. In this limit 
the wavelength is approximately given by A = iry/k, and again, the energy density 
is approximately constant. 

1.1 Domain wall solutions 

For k close to the critical value k c = 1 the solutions exhibit well separated wall- 
like structures. To illustrate this, we show in Fig.(|3]) a numerical solution <fi(x) 
of Class IA subject to the initial condition Re [0(0)] = 0.0001 and Im [0(0)] = 0. 
The value of the "constant of the motion" k that is associated with this solution is 
k 1 — 10~ 8 . The corresponding energy-density e(x) is also indicated. It is clear that 
if the trajectory of the solution in the complex plane passes the supersymmetric 
vacuum at the origin, the ratio of the width of the walls to distance between them 
is small since the force driving the trajectory away from the minimum is small. In 
Fig. (|) a numerical solution, belonging to Class II and subject to the initial condition 
Re [0(0)] = and Im [0(0)] = 0.0001, is shown. The value of the "constant of the 
motion" k that is associated with this solution is k « 1 + 10 -8 . For k j 1 the 
solution consists of equally spaced "right" walls (Class IA) or equally spaced "left" 
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1 2 3 4 5 

Figure 2: Wavelength A of the solutions as a function of k. 




Figure 3: Numerical solution to the BPS equation of Class IA subject to the initial 
condition Re [0(0)] = 0.0001 and Im [0(0)] = 0. Depicted are (a) the real part of 
4>(x), (b) the imaginary part of (f>(x), and (c) the energy density e(x). 

walls (Class IB), and for k j 1 the solution consists of equally spaced, alternating 
"left" and "right" walls. As the critical value k c — 1 is approached, the distance 
between the walls increases, but the width and the shape of the walls converge. 

1.2 Approximate solutions in between the walls 

In between the walls, both the real and the imaginary part of are very small. In 
this region the trajectory is approximately governed by the BPS equation linearized 
about = 

"202 , 

-20i . (9) 



dx 
dx 
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Figure 4: Numerical solution to the BPS equation of Class II subject to the initial 
condition Re [0(0)] = and Im [0(0)] = 0.0001. Depicted are (a) the real part of 
4>(x), (b) the imaginary part of 0(x), and (c) the energy density e(x). 

The general solution of Eq. @ is 

fa(x) = C+e 2x + C^e~ 2x , fa(x) = -C+e 2x + C.e~ 2x . (10) 

For Class I cycles, with the initial condition 0i(O) = u and 02(0) = 0, the solution 
is 

(f>i(x) = u cosh2x , 
02 (%) — — wsinh2x . (11) 

The constant u, which is very small for well separated walls, is related to the constant 
of motion k of the solution by k = (1 — u 2 /2) 2 . The approximation breaks down 
when x is so large that the magnitude of the fields 0i and 02 becomes of the order 
one. This happens at the location of the walls. The wavelength of the solution is 
therefore related to the constant k by A = — l/21n(l — k) for k | 1. 

Similarly, the approximate solution for Class II solutions in between a "left" wall 
and a "right" wall is given by 

0i (x) = — t>sinh2a;, 

4>2(x) = t>cosh2x. (12) 

The constant v, which also has to be very small for well separated walls, is related 
to the constant of motion by k = (1 + v 2 /2) 2 . From this approximate solution it 
follows that the wavelength of the Class II solution is related to the constant k by 
A = — ln(A; — 1) in the limit fell. 
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1.3 BPS domain walls on a circle 



So far we have discussed generic periodic solutions to the BPS equation. Now we 
want to determine what BPS saturated configurations are possible on a compact 
space with given circumference L. It is clear that we have to select only those 
solutions in the total solution space for which the wavelength A fits an integer N 
times on the circumference L. The allowed values of k are thus the solutions of the 
equation L = NX(k), where the integer N is a winding number that specifies how 
many times the solution winds around the poles when x varies between and L. In 
Fig.(0) we show the period function X(k), which gives the wavelength of a solution 
as a function of k. 

The simplest BPS configurations have winding number N — 1; for any value 
of L there is such a configuration of Class IA and Class IB, but only if L > X m 
is there such a configuration of Class II. The tension T = 2AW of the Class I 
configurations is equal to T = Att while the tension of the Class II configuration is 
equal to T = 87r. When L is very large, the Class IA and IB configuration describe 
just one "right" wall and one "left" wall, respectively. The Class II configuration 
contains one "right" wall and one "left" wall at a distance of L/2. 

Similar considerations apply to configurations with higher winding number N. 
The tension of the BPS configurations is proportional to N. As long as L is much 
larger than N times the width of the walls, the configurations contain well separated 
equidistant domain walls, N "right" ("left") walls in Class IA (Class IB) configura- 
tions, and N "left" walls alternating with N "right" walls in Class II configurations. 

1.4 Non-BPS domain walls 

As explained in the previous section there is no Class II BPS saturated configuration 
if the circumference L of the compact dimension is smaller than A m . However, there 
is a stable, static solution to the second order equations of motion with the same 
topology. Such a solution can not be BPS saturated. It is a non-contractable cycle, 
so it can not be shrunk to a point. At the same time it has to remain at finite 
distance of the origin in the complex plane because for fixed circumference L the 
kinetic energy grows if the length of the trajectory in the complex plane increases. 

Following the same argument, it is easy to see that there are stable, static con- 
figurations that are not BPS saturated in all other topological sectors. The simplest 
example is the homotopy class of non-contractable cycles which wind first clockwise 
around one pole and then counterclockwise around the other pole ("race-track" 
topology). With obvious notation, we will refer to these topologies as RL~ X and 
LR" 1 . ' 

Summarizing, in any non-trivial topological sector and for any value of L there 
is at least one static, stable solution to the second order equation of motion with 
minimal energy in that sector, a classical ground state. In certain topological sec- 
tors, those with the topologies R N , L N , (R _1 ) N and (L -1 )^, the stable solution 
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is BPS saturated for any value of L, and only half of supersymmetry is broken in 
the classical ground state. In other topological sectors, those with the topologies 
(RL) N and (i? -1 L _1 ) N , there are two stable solutions that are BPS saturated if the 
circumference L of the compact dimension is larger than a minimum value, NX m . 
If L is smaller than the minimum value, then there is still a static solution to the 
equation of motion, but it is not BPS saturated. 



2 d=3+l; localized light modes 

In this section we study the massless and light modes in a BPS saturated background 
that takes the form of a sequence of equidistant domain walls. We have discussed 
configurations of this type in Section [I]. As only half of supersymmetry is sponta- 
neously broken, the theory of the massless and light modes in d — 2 + 1 dimensions 
has M = 1 supersymmetry. We will therefore focus on the study of bosonic modes; 
clearly, each bosonic mode is accompanied by a fermionic mode that completes the 
supersymmetry multiplet. 

There is an exact bosonic zero mode originating from the invariance of the action 
under translations along the compact dimension of the BPS saturated background 
as a whole. We will now study the spectrum of the bosonic pseudo zero modes 
corresponding to independent translations of the domain wallsQ. 



2.1 Tight binding approach 

In oder to get a more quantitative picture of the bosonic pseudo zero modes, we 



follow Ref. |T3[ and expand the action in fluctuations of the fields around the clas- 
sical, BPS saturated background configuration up to quadratic order. In contrast 
to the case of a real wall background as it was investigated in Ref. |13[, we have 
to keep track of the fact that here the BPS saturated wall <p w is complex. Writing 



+ d>' and 



y — Yw T Y a ^ y — ^2 
of motion are obtained 



(d^ + WW + Re 
(d^ + W"W" - Re 



with real fields 0' x and 0' 2 , the following equations 



(13) 



WW" 


(f)'i — Im 


WW"' 


02 = 


WW" 


) 2 - Im 


WW" 


01 = 



The primes denote derivatives with respect to (0) in the case of W (W) evaluated 
at = <p w (0 = <j) w ). For localized massless and massive modes we make the 
following decomposition 



0i = Xi(z)A(t, x, y) , 02 = X2(z)A(t, x, y) 



(14) 



3 For large separation between the walls the action is nearly invariant under independent trans- 
lations of individual walls. 
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where x% an d A are real fields. From the equation of motion, Eg. (|13D , we obtain the 
following equations of motion for these fields 



8 n d n A {-d 2 + W"W + Re 



WW' 



Xi - Im 



WW 



X2 



A 

d n d n A (~d 2 z + WW - Re 



m 



Xi 
WW" 



X2 ~ Im 



WW 



Xi 



A 



X'i 



m 2 , (15) 



where n = 0,1,2 is a Lorentz index in d — 2 + 1 dimensions. Hence, the field 
A(t, x, y) has a mass m in d = 2 + 1 dimensions. The vector 



(16) 




X 

sat isfirs the eigenvalue equation 

where the hermitian operator H is defined as 



(17) 



H 



-efi + WW + Re 



— Im 



WW 



WW 



-Im 



WW 



-d 2 + WW - Re 



This equation has the form of a one-dimensional Schrodinger equation, and, in 
analogy, we will refer to H as the Hamiltonian and to x as the wave function. Using 
the BPS equation, Eq. (|7|), it is not difficult to explicitly show that H annihilates 
X = d z (p w , which corresponds to the Goldstone boson of spontaneously broken overall 
translational invariance along the z direction. 

Since we are not able to solve the eigenvalue problem in Eq.([L7]) exactly, we resort 
to the tight binding approximation to study the masses of the pseudo zero modes 
corresponding to independent translations of the individual walls. For definiteness, 
we take the R topological sector, with the background of iV walls separated by a 
distance A = L/N. Considerations for the (RL) N sectors are similar, but slightly 
more complicated due to the sequence of alternating walls and anti-walls. We will 
uncover the light bosonic modes of the wall lattice by first considering a massless 
mode bound to a single wall in isolation, just as one can find the electronic levels of 
a one dimensional crystal from the energy levels of the individual atoms when they 
are far apart and overlaps are small. 

Thus, we consider a single wall, </>^, along an infinite dimension. In this limit 
is a normalized eigen function with m 2 = 0. In the case of N such 



X 



walls separated by an infinite distance, there are iV localized (bound) wave functions 
with vanishing mass. However, when the walls are at a large, but finite, distance 
from each other, the iV-fold degeneracy is lifted due to the small overlap of potentials 
and wave functions centered at different wall sites, and a band structure emerges. 
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We will now elucidate this band structure, and we will estimate the width of the 
band. 



To this end we decompose the operator H c = H° + V c , where 



-&l + 4 
-dl + A 



and 



V c 



W"W" - 4 + Re 



WW 



-Im 



ww: 



— Im 

W'JW'J - 4 



ww 

-Re 



WW 



(19) 



(20) 



Here the superscript c indicates that the quantity is evaluated at <p = <f>™. With this 
definition, V c vanishes exponentially fast away from the wall. We approximate H 
in the case of N walls separated by a distance A as 



(21) 



This Hamiltonian commutes with the generator of the Rjy symmetry, the transla- 
tions over distances nX. The (pseudo zero mode) wave functions, the equivalent 
of Bloch waves, are linear combinations of x° centered at each of the N wall sites. 
These linear combinations transform under irreducible representations of Rn- They 
take the form 

Xk = £a fc [n]x°(*-nA), (22) 



where a&, with k 
matrix 



1..N, are the N real, orthonormal eigenvectors of the N x N 



( 2 



-1 
2 -1 
-1 2 







1 \ 



-1 
















V 











(23) 



2 J 



This matrix has a non-degenerate eigenvalue equal to zero for any value of N, 
with eigenvector d\ = t^(1> 1, •-, 1). For even values of iV there is a second non- 
degenerate eigenvalue equal to 4 with eigen vector a 2 = 7^(1> — 1, 1, •-, — 1, 1). All 
other eigenvalues are between and 4 and are doubly degenerate. Therefore, the iV 
dimensional representation of R^ is reduced to 2 one- dimensional and N/2 — 1 two- 
dimensional irreducible representations in case iV is even, and to 1 one-dimensional 
and (N — l)/2 two-dimensional irreducible representations in case N is odd. 

As an aside, the eigenvectors and eigenvalues of the matrix in Eq. (|23|) also yield 
the normal modes and frequencies of a closed linear chain of iV identical masses and 
springs. In d = 3 + 1 dimensions, the walls extend into two infinite dimensions and 
have therefore infinite mass. However, for solitons in d = 1 + 1 dimensions, which 
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have finite energy, these vibrations are physical. If the distance between the solitons 
is large as compared to their width, these vibrations are very soft. 

Given the wave functions in Eq. fl2"2]) , the approximate mass eigenvalues are then 
given by 

= It H* 

JO A 

In fact, for the wave functions that transform under the two dimensional representa- 
tion of Rn, degenerate perturbation theory should be used. Including only overlap 
integrals from neighboring sites, and using the asymptotic form for x° far away from 
the wall, we estimate that the masses are exponentially suppressed, 

m 2 TB = Ce- 2X , (25) 

where C is of the order one. We thus see that the N degenerate vanishing masses 
in the limit where the walls are infinitely far apart evolve into a (mostly double 
degenerate) band of masses when the distance between the walls becomes finite. 



2.2 Oscillations 

For an observer who can not resolve the circumference L of the compact dimension, 
the physics of the light modes is described by an M = 1 supersymmetric theory in 
d — 2+1 dimensions with a spectrum as we discussed in the previous section. In such 
a world there is a hierarchy between the heavy modes, with masses comparable to the 
wall tension, and the light modes, with masses that are exponentially suppressed by 
the ratio of the distance between the walls and their width. It is interesting to note 
that in topological sectors with multi-wall backgrounds that are not BPS saturated 
the physics of the light modes is described by a M = 1 model with supersymmetry 
breaking. The scale of the breaking is equal to the scale of the light masses. It 
seems that under such conditions there are no hierarchy or naturalness problems. 

Alternatively, an observer, who is confined to a particular wall and can not resolve 
its width, experiences a different physical reality. Such an observer does not have 
the means to do experiments in which the massive bulk modes are excited. Still, due 
to the BPS saturation of the background, his world has M = 1 supersymmetry. But 
he will be perplexed by the appearance and disappearance of particles in his world. 
From the ansatz of Eq. ([22] ) it is clear that a state localized on a particular wall must 
be a superposition of the wave functions that represent the light modes. This in 
turn means that such a state is not a mass eigenstate. The implication is oscillation 
of localized particles between the various walls with frequencies corresponding to 
the mass differences. In the analogy of a linear crystal, this would correspond to 
electrons hopping from atom to atom. For example, in the case of the R 2 sector 
with two walls, the approximate massless mode is 

&{z)=j?(z)+j?{z-\) , (26) 
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and there is one light mode, approximately given by 

Uz) = *°(*) - f(z - A) . (27) 

A state localized one of the walls is constructed by the superposition of the mass 
eigenstates, 

Xo(z)±Xi(z)- (28) 

3 d=2+l; chiral fermions 

When the model is dimensionally reduced to d — 2 + 1 dimensions, it has M = 2 
super symmetry. It has inherited four supersymmetry generators, whereas minimal 
supersymmetry in d — 2 + 1 requires only two generators. The scalar potential may 
be written as 

-4M 

The object W plays the role of the superpotential in the dimensionally reduced 
model. It is related to the superpotential W in d = 3 + 1 dimensions by 



W = 21m 



^(t2 (01 + ' 02) 



2 arctan — ^ 2 = + 2 arctan — ^ 2 ■= . (30) 
01 + v^ <j>i-V2 



The extended supersymmetry manifests itself in the fact that the superpotential W 
is harmonic. In terms of the scalar components of the Af = 1 superfields in d — 2 + 1 
dimensions, 0i and 2 , the BPS equations are 

dth, d\V 

a— — = 2a 



dx d(j>i " \ (0 X + y/2) 



2 + 0i 




dx a d<p 2 a \(0 1 + V2) 2 + 0| (0 X - V2y + <P 2 J ' 1 } 

where a — ±1. The interesting new feature in d — 2 + 1 dimensions is that now 
"chiral" fermions exist on the d — 1 + 1 dimensional worlds of the domain walls. In 
d — 1 + 1 dimensions there is no spin, but massless fermions can be either left or right 
moving; this is what we will allude to as "chirality" . As outlined in the Appendix, in 
a BPS saturated background configuration there is always one left moving and one 
right moving fermionic zero mode, corresponding to the two spontaneously broken 
supersymmetry generators. Whether a mode that is associated with a particular 
broken generator is left or right moving depends on the sign a in the BPS equation 
that provides the background configuration. In other words, it is not important 
around which pole the background configuration winds in the complex plane, but 
whether it winds clockwise or counter clockwise. 
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3.1 Wess-Zumino model 



In order to further investigate the issue of the "chirality" of massless fermions on 
the wall, we discuss in this section the regular Wess-Zumino model reduced to d — 
2 + 1 dimensions. It is well known that this model allows for real, BPS saturated 
domain walls (kink or anti-kink). However, these solutions do not exist in a compact 
dimension, and there are no BPS saturated multi-wall solutions (although almost 
static configurations with walls at large distances do exist). For the issues we want 
to discuss here, these differences are irrelevant. 

Our analysis is similar to the discussion of fermionic zero modes in Ref . Hl"9| , 



where the same model was investigated in d = 1 + 1 dimensions. The superpotential 
of the dimensionally reduced, regular Wess-Zumino model is 

W = 72 (^1-^ + ^2)- (32) 

There are two vacua, 0i = ±\/2 with <p 2 = 0. The BPS saturated domain walls 
connecting these vacua are determined by 



dy v/2 

^ = 0^(20^), (33) 

where a = ±1. The solutions that interpolate between the two vacua (the well- 
known kink or anti-kink, depending on the sign a), take the form 

0™ = ocy2 tanh(x — Xq) 

02 = 0. (34) 
The equation of motion for the fermions in the background is 

il m d m i) t - W l3 ^ = , (35) 
where the "mass" matrix W is given as 

There are two massless fermions on the wall corresponding to Rebbi-Jackiw zero 
modes, 

ijji{x,y,t) = A(y)((x,t), 

il)*{x,y,t) = 0, (37) 



with 



A(y) = e J v . 



ry 2 £ = a( 
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and 



0, 

B(y)i](x,t) , 



(38) 



where 




e 



—a 



£ w 2 2(y'W 

J yn 



—ar] . 



Both ( and r\ satisfy the massless Dirac equation in d = 1 + 1 dimensions, ry a <9 a £ = 0, 
where a = 0, 1. Note that 7 2 is the chiral projection matrix in d = 1 + 1 dimensions, 
so that one of the massless fermions is right moving and the other one is left moving. 
These massless fermions are the goldstino fermions associated with the breaking 
of two out of four generators of the Af = 2 supersymmetry, as discussed in the 
Appendix. 

We proceed by addressing the fate of these massless fermions when the Af = 2 
supersymmetry is explicitly broken to Af = 1. In this case, the wall breaks one 
of the remaining two supersymmetry generators, and therefore only one Goldstone 
fermion exists. What happens to the second massless fermion when the explicit 
supersymmetry breaking is switched on? To answer this question we introduce 
explicit Af = 2 breaking terms in the superpotential. 

First, we insert a coupling constant A in front of the 0i0 2 term in the superpo- 
tential, so that the Af = 2 supersymmetry is explicitly broken unless A = 1. The 
vacuum expectation values and the wall solution do not depend on A, but W 2 2 now 
is equal to X(j)f. Even though there is only one Goldstino, from the perspective of 
the Rebbi-Jackiw zero modes in Eqs. ( f3?l , |38|) , it is clear that the second massless 
fermion continues to exist, even if A ^ 1. When A changes sign, the second massless 
fermion changes from right to left moving (or vice versa, depending on the sign a), 
and for A = the second massless fermionic mode is not normalizable. 

Next, we add a term | to the superpotential, so that the Af = 2 supersym- 
metry is explicitly broken to A/" = 1 for fi ^ 0. Again, the vacuum expectation 
values and the wall do not depend on /i, but now W 22 = <pf + In this case the 
second massless fermion persists as long as < \/2. For > \/2 the second zero 
mode becomes non-normalizable. 

It is clear that since the fermionic modes on the wall are chiral, they can only 
become massive in pairs, combining one left moving and with one right moving 
fermion. Individual massless modes, either left or right moving, can only disappear 
when they become non-normalizable. 

3.2 Complex, mult i- wall backgrounds 

In the case of a BPS saturated background with one domain wall along a compact 
dimension, there is one massless scalar localized on the wall corresponding to the 



16 



spontaneously broken translational symmetry. Moreover, there is one left moving 
and one right moving massless fermion, corresponding to the two spontaneously bro- 
ken supersymmetry generators (see the Appendix). These massless modes combine 
to form a massless superfield in d — 1 + 1 dimensions. 

The supersymmetry multiplet in d — 1 + 1 dimensions contains a real scalar 
and a Majorana fermion. Therefore, in the case of N walls, there will be a total 
of iV — 1 light superfields (d = 1 + 1, M — 1 supersymmetry) with exponentially 
suppressed masses, and one massless superfield. In order to make states localized 
on a particular wall, states with different masses have to be superimposed, just as 
in the d = 3 + 1 dimensional case. 

The explicit analysis of the Rebbi-Jackiw modes becomes impossible if the 
fermion "mass" matrix in the background is not diagonal, as is the case with the 
type of models we consider in this paper. But the results of the previous section 
must still be valid. It is therefore expected that if the M = 2 supersymmetry is 
explicitly broken to M = 1, there will still be a left and a right moving fermion 
on the wall, at least for a finite range of the breaking parameters. Moreover, in a 
compact dimension, it seems unlikely that fermionic zero modes can disappear by 
becoming non-normalizable. 

4 d=l+l; solitons 

The main new feature in d = 1 + 1 dimensions is that the domain walls, which have 
infinite energy in higher dimensions, now reduce to solitons with a finite mass. Such 
solitons have a particle interpretation, whereas the domain wall backgrounds in the 
higher dimensions are interpreted as vacua. 

All topological sectors in the d = 1 + 1 dimensional theory, whether they have 
a BPS saturated ground state or not, can be assigned a pair of winding numbers 
(Ni,N r ), where iVj e Z is the winding number around the "left" pole at <fi = —1 
in Fig.(|IJ), and AV G Z is the winding number around the "right" pole at <fi = 
+ 1. According to our conventions the winding number increases by one if a pole 
is encircled in the counter clockwise direction. This classification of the sectors 
according to the pair of winding numbers is not unique; for example, the sectors 
RLRL' 1 and R 2 have the same winding numbers, (0,2), but they are topologically 
distinct. 

The sectors L 1 and R T have a classical ground state in which only half of 
supersymmetry is broken for any value of L. In the L 1 sector this ground state 
corresponds to the solution of the BPS equation which in the large L limit contains 
Ni "left" solitons if Ni is positive, or |JVj| anti "left" solitons if Ni is negative, 
separated by a distance L/\Ni\. The situation is analogous for the R Nr sectors. 

For fixed L > A m , the sectors (RL) N have two degenerate classical ground states 
in which only half of supersymmetry is broken. In the large L limit, and if N is 
positive, one of these ground states corresponds to the solution of the BPS equation 
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containing N "left" solitons which alternate with N "right" solitons, separated by 
a distance L/2N. Similarly, if N is negative we have the BPS solution with \N\ 
"anti-left" solitons alternating with \N\ "anti-right" solitons separated by a distance 
L/2N. The second, degenerated ground state corresponds to the BPS solution 
that winds around both poles at large distance, iV times in the counterclockwise 
direction if N is positive, and \N\ times in the clockwise direction if N is negative. 
This solution to the BPS equation has approximately constant energy density in the 
large L limit. There is a tunneling transition with finite action that connects the two 
ground states and the remaining 1/2 of supersymmetry is broken non-perturbatively. 
For L < NX m , there is no BPS saturated ground state in the sector (RL) N already 
at the classical level. There are static, stable solutions to the second order equation 
of motion which break supersymmetry completely at the classical level. 

For the same reason, supersymmetry is broken completely in the groundstate 
at the classical level in all sectors other than the ones we discussed above. In all 
sectors, the ground state energy tends to 2ti(\Ni\ + |iV r |) in the large L limit. The 
theory has two types of solitons, "left" and "right" , and two topologically conserved 
quantum numbers, Ni and N r . The solitons interact through a short range force, 
which decreases exponentially with distance. "Left" and "right" solitons repel each 
other, and so do "anti-left" and "anti-right" solitons. "Left" and "anti-left" solitons 
on the other hand attract each other. 



4.1 Semi-classical calculation of the tunneling amplitude 

In this section we calculate the amount by which the ground state energy is lifted 
above the BPS bound by the tunneling transition in the sector (RL) N , for N = 1 
and L > A m . We choose N = 1 for definiteness; the tunneling action scales linearly 
with \N\, and therefore the result can be trivially generalized for any value of N. 
The two ground states correspond to two solutions of the BPS equation in Class II, 
one with k < k m and one with k > k m , as shown in Fig.([5|). At the non-perturbative 
level the corresponding solitons are not longer BPS saturated because a tunneling 
process mixes them and lifts their mass above the BPS bound. This phenomenon was 



first observed in Ref . ||20|| . In Ref.|2T| it was noted that the calculation of the semi- 
classical instanton action is equivalent to the calculation of the energy of domain 
wall junctions in d = 2 + 1 dimensions. In the same work, the general calculational 
framework of the action was established and illustrated for a specific model. 

We assume that the instanton in the imaginary time saturates the BPS bound 
in the model we consider; there is no obvious reason why this would not be the case. 
If the instanton is BPS saturated, it satisfies the equation 

^ - i e «^ (39) 

where £ = x + it, and t is the euclidean time. The field <fi depends both on ( and (. 
The phase 5 has to be chosen so that it is consistent with the boundary conditions 
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Im[ 4f] 




Re[ 4f] 



Figure 5: The initial (outer) and final (inner) configurations in the complex <fi plane. 
Both configurations are solutions to the BPS equation with the same wavelength 
A. The instanton interpolates smoothly between the initial and final configuration 
with finite action. 

that are imposed on the solution (see below). The tunneling action takes the form 

+ / dxdt—-— T . (40) 



In Figs. (|^, ||) we have indicated the boundary conditions that are imposed on the 
instanton configuration. At t — —T/2 the instanton field <ft(x,t) takes the form 
of the "outer" configuration (see Fig. (j^) ) and at t — T/2 it takes the form of 
the "inner" configuration. We always have in mind the limit T — > oo. Periodic 
boundary conditions apply to the vertical edges in Fig.(|6|). We have chosen the 
solutions to wind counter-clockwise in the <p\-(f)2 plane when followed in the positive 
x direction, which corresponds to 5 = 7r/2 in the BPS equation. The initial and final 
configuration are positioned such that <f>(L/2, — T/A) = ivo and 0(L/2,T/4) = ivi. 
Note that in order for the instanton configuration to satisfy the BPS equation, it has 
to connect the solution with the larger value of k as the initial configuration with the 
solution with the smaller value of k as the final configuration. The anti-instanton 
that interpolates in the reverse situation satisfies an anti-BPS equation. The action 
can now be written as 



Abps = / dxdt 



2V'5-(Vxa) 2 , (41) 
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Figure 6: Boundary conditions in the x-t plane. The initial and final configuration 
at the bottom (—T/2) and top (T/2) of the rectangle are traced counter-clockwise 
in the complex plane (see Fig.(|5|)). Periodic boundary conditions apply to the 
vertical edges. The shading schematically indicates the energy density, with the 
darker areas indicating higher energy density. 



20 



where 



and 



cos 5 Re [W] + sin<5Im [W] 
costflm [W] - sin 5 Re [W] 



z / wsu rus [vv j -r sin u uii [vv j \ //io\ 



(43) 



With the application of Gauss' and Stoke's theorems the surface integrals can be 
converted into line integrals over the boundaries of the surface, 



/ Im 


'*% 




\ Im 




) 



A 



BPS 



2 j> S ■ dn 



a ■ dx, 



(44) 



where the boundaries are the edges of the rectangle in Fig.(|B|). Due to the periodic 
boundary conditions only the integrals of S contribute over the vertical edges of 
the contour, as the superpotential is multi-valued. Turning now to the explicit 
calculation of the tunneling action, there are three non-vanishing contributions. 
First, consider 



4° 

A BPS 



l-T/2 

2| r/2 {lm[^L =# -Im[iyL = _ # }rft, 



8ttT. 



(45) 



This contribution is equal to the classical ground state energy times the elapsed time 
T. It is equal to the action in case there is no tunneling, when the initial and final 
configuration are the same. We define the tunneling probability by subtracting this 
trivial contribution in the exponent. Next, consider 



BPS 




{Re[W] t= _T-Re[W] t= r}dx 



(46) 



where k Q and ki are the two solutions to the equation X(k) = L with k > k c . We 
define k a to be the larger of the two. Note that the integrands are "constants of 
motion" for the solutions to the BPS equation, and therefore the integrals can be 
expressed in k Q and ki, respectively. In the limit A = L ^> A m , ki approaches 1 from 
above, and k Q is related to L by L = ny/k. In this limit, the second contribution to 
the tunneling action is therefore A l BPS = 2LIiiL/tt. Finally, we calculate 



A 2 

^BPS 



1 j Im 




— Im 




J 




dx 


1 2 


dx 





dx 



-S{k ) + S(h 



(47) 



where S(k) is the area enclosed in the 0i-02 plane by the solution to the BPS 
equation marked by k. S(k) is a monotonically increasing function of k. In the limit 
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k I k c = 1, S(k) approaches S(l) = 4 from above. For large values of k on the other 
hand S(k) = 2-K\fk. Therefore, in the limit A = L ^> A m , the third contribution 
to the tunneling action is approximately A 2 BPS = —2ny/k^ + 1, or, in terms of L, 
A BPS = -2L + l. 

We have calculated the tunneling action solely in terms of quantities related to 
the boundary configurations. An explicit interpolating instanton configuration is not 
needed, just the existence of a BPS saturated instanton configuration is sufficient. 
In the large L limit, the semi-classical approximation is valid, and the ground state 
energy is lifted above the BPS bound by an amount 

AE oc 8ne 2ABPS = 8vr e~ 4LlnL , (48) 

to leading order in L. The solution to the BPS equation that connects the vacuum 
at <fi = to the runaway vacuum at — > ioc (this is a solution that is not periodic, 
obtained with 5 = 0) has infinite tension. The circumference L in a sense acts 
like an infra-red regulator. For this reason, the tunneling action is not linear in L, 
but proportional to LlnL. We also note that in this model the instanton mixes a 
state containing two localized solitons with a state that corresponds to a flat energy 
density. This is in contrast with the model considered in Ref.|21]], where the mixing 
is between two states with different solitons. 

Summary 

We have constructed a toy model to study some issues related to the idea of "brane 
worlds" in a field theoretical context. This toy model is a simple Wess-Zumino 
model with only one chiral superfield in d = 3 + 1 dimensions. Different aspects 
were studied in various dimensions. 

We showed that it is possible to have BPS saturated configurations with multiple, 
well separated domain walls in a compact space dimension. Generically, in such a 
background there are light modes with masses exponentially suppressed by the ratio 
of the distance between the walls and their width. These modes form the "matter" 
of the "brane world" , wheras the domain walls play the role of branes in our simple 
picture. We also described how particles can oscillate between walls. We illustrated 
these issues in our specific model. 

We also found that the "chiral" fermions localized on d = 1 + 1 dimensional BPS 
saturated walls embedded in a d = 2 + 1 dimensional space-time always come in 
pairs; one is left moving, and one is right moving. 

Finally, we discussed the toy model in d = 1 + 1 dimensions. The domain walls 
then reduce to solitons with a finite mass. We showed that some classically BPS 
saturated states are lifted non-perturbatively by instanton tunneling. 
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Appendix A 

In this Appendix we establish our notation and discuss the zero modes in a BPS 
saturated background. 

d=3+l dimensions 



We use the conventions outlined in | 13| . The signature of the metric g^ u is (+, — , — , — ) 
The component Lagrangian for an M = 1 Wess-Zumino model in d = 3 + 1 dimen- 
sions with canonical Kahler potential and superpotential W is 

dW ld 2 W in , 

-dW ld 2 W - T . , , 

+F W " 2~&f (49) 

The component fields transform under the supersymmetry transformations as 

5F = i(a^ &a d^ a , (50) 

where the transformation parameter ( is the two component complex Grassmann 
variable._ The auxiliary field F can be eliminated by the use of its equation of motion 
= — ^ to give the Lagrangian density 

ld 2 W ld 2 W - -■ 

-YW^-TW^'"- (51) 

The BPS equation for domain walls reads 
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where in a non-compact direction the phase 5 is determined from the equation 



Im 



e- td W(<p(z = -oo) = Im e~ l8 W{<\){z = +00) . (53) 



BPS domain walls in a compact dimension are only possible if the superpotential is 
not single valued. In the models we consider e l& ± i. 

Let us now expand the Lagrangian of the model about the non-trivial solution 
4> w of the BPS equation. Setting = <p w + 0' and expanding up to quadratic terms 
in the fluctuations 0', we obtain 

{w"'W'(p' 2 + W'W'y' 2 + 2W"W"(f)'4>' + W"ip a ip a + W'i)^} ,(54) 

where the primes on W (W) indicate derivatives with respect to (0) evaluated at 
= (fiw (0 = <t> w )- The equation of motion for the scalar in the wall background is 

<9 2 0' + {W'W'"$ + W"W"(j)'} = 0, (55) 

and the equation of motion for the fermion in the wall background is 

<Ar - W"ij a = 0. (56) 

Half of the supersymmetry generators is broken by the wall solution. In order to 
determine which generators are unbroken, we act with the SUSY generators on the 
wall and find 

^« = ~^(C + ««^C 4 ), (57) 

where we have used the BPS equation, and a — ±1. Transformations, for which the 
parameter satisfies ( a = i' a(J aaC a i correspond to unbroken symmetries. The massless 
fermions localized on the wall which are related to the broken SUSY generators take 
the form 

ip a (x, y, z, t) = --^f( z )(c<(x, V, t), (58) 
where ( a (x,y,t) is constrained by 

C a = +taa^C- (59) 

It can then be verified that ip a satisfies the fermionic equation of motion in the 
domain wall background if ( a satisfies the equation 

t^ a a^d m ( p = 0, (60) 

where m takes the values 0,1,2. This is exactly the Dirac equation in d — 2 + 1 
dimensions, 

l^dmC = 0, (61) 
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with the following representation of the 7 matrices; 7 = a 3 , 7 1 = — 1a 2 and 7 2 = io\. 
In this basis the charge conjugation matrix, defined by C T = —C and C^^C^ 1 = 
— (7 At ) T , takes the form C = ±ta 2 . The Majorana constraint reads 

( T C = (62) 

This is identical to the constraint in Eq.(^). We thus find a massless Majorana 
fermion in the d = 2 + 1 dimensional world on the wall. Together with the massless 
real scalar, this fermion forms a SUSY multiplet in the lower dimensional theory. 



d=2+l dimensions 

M = 1 SUSY 

The Lagrangian of the model with M = 1 supersymmetry and one superfield with 
canonical kinetic terms and superpotential W takes the form 



£ = \{d m <pd m <P + i^ m d m ^ + F 2 



<9W d 2 W - ) 

+2F ^-^\- < 63) 

We choose the Majorana basis for the 7 matrices, 7 = a 2 , 7 1 = 1a 3 and 7 2 = 1a 1 . 
Here is a real scalar field and ip is a rea l two component spinor, with the Dirac 
conjugate defined by ip = ip Tr y° . The supersymmetry transformations take the form 

5^ = -*<9 m 7 m C + , (64) 
5F = -iCl m d m ^ , 

where the transformation parameter ( is a two component real Grassmann variable. 
The BPS equation is 

dy d(p 

where a — ±1. 

Expanding the Lagrangian about the non-trivial solution <p w to Eq. fl65"D up to 
quadratic terms in the fluctuations, we obtain 

£ = 1 [d m (j)'d m 4>' + z^y"^ - (W" 2 + W'W'") 0' 2 - W'lpip} . (66) 
The equations of motion for the scalar and the fermion are 

d 2 <f)' + {W 2 + WW") 0' = 0, (67) 

and 

i~1 m d m i) - W"ip = 0. (68) 
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The scalar zero mode related to the broken translational symmetry in the y 
direction is 

tf/(x,y,t) = ^-(y)A(x,t), (69) 

where A(x,t) satisfies the equation of motion of a massless real scalar, d 2 A = 0, in 
d — 1 + 1 dimensions. By acting with a supersymmetry generator on (fi w we find 

<ty = - a *pL(i - m7 2) C . (ro) 
dy 

For the unbroken generator ( = 1/2(1 + cw7 2 )C, whereas for the broken generator 
( = 1/2(1 — mj 2 )(. The fermionic zero mode related to the broken supersymmetry 
generator is 

Mx,y,t) = y±(y)ax,t), (71) 

where ( satisfies the massless Dirac equation in d — 1 + 1 dimensions ry a <9 a C = 0, 
with a = 0, 1, satisfies the Majorana constraint and in addition the constraint ( = 
1/2(1 — cn7 2 )C, that is, ( is chiral. Consequently, the massless fermion on the wall 
is either left or right moving depending on the sign of a. 

N = 2 SUSY 

The minimal model with M = 2 supersymmetry in d — 2 + 1 dimensions contains 
two Af = 1 supermultiplets. The Lagrangian in components is 



C = \{d m <p i d m <p i + iA 1 m d m A + 



dW <9 2 W - 1 

where the indices i, j are 1, 2 and the superpotential is harmonic, so that 

d 2 W d 2 W 



<90? d<p 2 



0. (73) 



An J\f = 2 model in d — 2 + 1 dimensions can be obtained from an A/" = 1 model 
with one chiral superfield in d — 3 + 1 dimensions. In this procedure, the fields are 
assumed to be independent of one of the spatial coordinates. The superpotential 
in the lower dimensional theory is obtained from the superpotential in the higher 
dimensional theory by 



W(0i,0 2 ) = 2Re 



(74) 



The lower dimensional superpotential obtained in this way is automatically har- 
monic, because the higher dimensional superpotential is holomorphic. The N = 2 
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model in d = 2 + 1 dimensions has, apart from the regular supersymmetry transfor- 
mations of N = 1 models, 

<fyi = -^ m l7 m C + ^iC, ^ = -«^27 m C + ^C, (75) 
5Fi = -<7 m <9«^i , SF 2 = -iCl m d m ip 2l 

additional supersymmetry transformations that take the form: 

<50i = -77-02 , ^02 = #1 , 

5-01 = -td^j^r] - F 2 r) , = ^017^ + ^i?7 > (76) 

From the perspective of dimensional reduction, these extra transformations arise 
because the underlying d = 3 + 1 dimensional theory has four supersymmetry gen- 
erators. 

The auxiliary fields Fi can be eliminated to give the Lagrangian 
The BPS equation is 

#i = a— (78) 
rfy <90j ' 

with a = ±1. 

Expanding the Lagrangian of Eq. ( |77| ) about the solution 0™ up to quadratic 
fluctuations, we obtain 

£ = \ {d m <P' t d m <P' t + i$a m d m ^ - WiWmM'i - WikWu&'A - >V,, r, r., } , (79) 
where the subscripts on W indicate derivatives with respect to 0j evaluated at 

& = 0r- ' 

In the domain wall background, the scalar equation of motion is 

d m d m <p' n + {WiW M + WinWu} 0; = 0, (80) 

and fermions satisfy 

i^d^ - Wijiftj = 0. (81) 

There is a massless scalar, A(x,t), on the wall which is related to the broken trans- 
lational symmetry, 

dd> w 

ti(x,y,t) = -£-(y)A(x,t) , (82) 
dy 

and a massless fermion, ((x,t), (either right or left moving, depending on the sign 
of a), due to the broken regular supersymmetry generator, 

d6 w 

V^-^COM), (83) 
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where ( is a (real) two component Major ana spinor satisfying the 1 + 1 dimensional 
massless Dirac equation and the constraint ( = 1/2(1 — ia r y 2 )C > . In addition, one of 
the extra SUSY generators is broken, which results in another massless fermion on 
the wall, 



ipi{x,y,t) = 




Mwt) = --^rj{x,t), (84) 

where rj is another (real) two component Majorana spinor satisfying the 1 + 1 di- 
mensional Dirac equation and the constraint ( = 1/2(1 + ia , y 2 )(. We therefore find 
in the case of the BPS wall in the d = 2 + l, J\f = 2 massless scalar, a left 

moving massless fermion and a right moving massless fermion, which together form 
an M = 1 multiplet in the lower dimensional theory. 
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